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The recent high statistics measurement of thepp — > ppr] reaction at an excess energy Q=15.5 MeV 
has been analysed by means of partial wave decomposition of the cross section. Guided by the dom- 
inance of the final state ^So pp interaction (FSI), we keep only terms involving the FSI enhancement 
factor. The measured pp and r;p effective mass spectra can be well reproduced by lifting the standard 
on-shell approximation in the enhancement factor and by allowing for a linear energy dependence in 
the leading ^Po — > ^So, s partial wave amplitude. Higher partial waves seem to play only a marginal 
role. 
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I. INTRODUCTION 

In recent years major advances have been made in 
the experimental 0, 0, IS S IS and theoretical 
IS IS kPj 10,, iSl investigation of the near threshold me- 
son production reactions in nucleon-nucleon collisions 
(for a comprehensive review cf. 111). With the ad- 
vent of medium-energy accelerators (ICUF, CELSIUS 
and COSY) and the corresponding influx of high pre- 
cision data on the total and differential cross section as 
well as the polarization observables, it has become possi- 
ble to study the interaction of the flavor-neutral mesons 
(eg. tt'^, and 77') with nucleons. A direct insight into 
such interaction can not be gained from meson-nucleon 
scattering experiments as the latter are impractical ow- 
ing to the very short life time of these mesons. Naturally, 
the largest data base has been accumulated for the pions 
but the bulk of data on 77 production inproton-proton 
collisions has also expanded significantly |S IS HI - The 
rj meson, which is the next lightest non-strange member 
of the pseudoscalar octet, has focused considerable at- 
tention of the nuclear community since it was established 
that the 77-N interaction was quite strong and attractive 
which might lead to a possible existence of the 77-nuclear 
bound states. 

In the recent measurements of the pp — > pprj reaction 
a very accurate determination of the four-momenta of 
both outgoing protons allowed for the full reconstruction 
of the kinematics of the final rjpp state. In consequence, 
these measurements provided in addition to the 77 and 
the proton angular distributions, also the pp and rjp ef- 
fective mass distributions 3, IS ■ The common feature 
of the near-threshold meson production in proton-proton 
collisions is the dominance of the very strong proton- 
proton final state interaction. The Monte Carlo simula- 
tions as well as direct calculations reveal that for small 
excess energies the inclusion of FSI enhances the cross 
section by more than an order of magnitude. This effect 
is also clearly visible in the effective mass distributions: 
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as a prominent peak close to threshold in the pp effec- 
tive mass distribution or as a bump near the end-point 
in the rjp effective mass distribution. The description in 
terms of a simple model in which a constant rj production 
amplitude is multiplied by an on-shell proton-proton FSI 
enhancement factor is only qualitatively correct, but at 
a quantitative level turns out to be insufficient to explain 
the experimental pp and Typ effective mass distributions. 
The contribution from higher partial waves and the final 
state 77P interaction have been indicated as the possible 
sources of this discrepancy. In comparison with pp inter- 
action, however, the role of 77-N interaction is much less 
important and taking as a rough measure the ratio of the 
77P and pp scattering lengths squared we may expect an 
effect at the level of about 1%. A full three-body calcula- 
tion based on the hyperspherical harmonics method ac- 
counting for all pair-wise final state interactions corrobo- 
rates the above estimate showing that the 7/p interaction 
modifies significantly only the total 77 production cross 
section as function of the excitation energy whereas the 
distortion of the shapes of the effective mass distributions 
does not exceed the current experimental uncertainties. 
The hyperspherical harmonics method will be presented 
elsewhere [TJ and here instead we would like to confine 
our attention to a purely phenomenological description 
of the effective mass spectra at the lowest available exci- 
tation energy equal Q=15.5 MeV. At such a low Q only 
a few partial waves in both initial and final states are ex- 
pected to participate in the 77 creation process which sub- 
stantially facilitates the interpretation of the experimen- 
tal findings. Among the possible amplitudes the partial 
wave amplitude in which the two protons are deposited in 
the ^5'o state plays a dominant role because it is propor- 
tional to the large FSI enhancement factor. Therefore, in 
the first approximation, in the cross section it should be 
sufficient to retain the square of the modulus of the dom- 
inant amplitude plus the interference terms that would 
contain the large FSI enhancement factor. However, after 
integrating over the angles the interference terms vanish 
and give no contribution to the pp effective mass dis- 
tribution. This implies that there is really no room for 
improvements left unless a weak energy dependence is 
admitted in the dominant partial wave amplitude. We 
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demonstrate in this paper that with this rather modest 
assumption both pp and rjp effective mass distributions 
can be quite weil reproduced. 

The pian of our presentation is as follows. In the next 
Section the FSI enhancement factor is revisited. We ar- 
gue that owing to a too steep a fall of the on-shcU en- 
hancement factor the approximate form thereof should be 
abandoned if favor of the full off-shell expression. Hav- 
ing established the best form of the enhancement factor 
we discuss the partial wave expansion of the pp — > rjpp 
transition amplitude in a quest of an approximate ex- 
pression for the cross section that would be valid for the 
lowest excitation energies. Finally, in Section III we ver- 
ify our simple model by presenting the comparison with 
experiment. 



II. THEORETICAL FRAMEWORK 

A. Derivation of the enhancement factor 

Since the proton-proton final-state interaction is be- 
lieved to be the dominant ingredient in the description of 
the pp pprj reaction close to threshold, it is logical that 
we first wish to reexamine the FSI problem to make sure 
that adequate measures have been taken to obtain the ul- 
timate solution. The basic idea how to account for final 
state interaction was put forward 50 years ago by Fermi 
1 Watson [T^ , Migdal |T5| and others (for a review cf. 

is based on the observation that in many 
processes the interaction responsible for carrying the sys- 
tem from the initial to the final state is of such a short 
range that in the first approximation may be regarded 
as point like. As a prototype one may consider a meson 
(x) production reaction NN NNx. To generate the 
meson mass m in nucleon-nucleon collision a large mo- 
mentum transfer is required between the initial and the 
final nucleons, which is typically of the order vMm, with 
M being the nucleon mass. The corresponding "range" 
of the production interaction is therefore much shorter 
than the range of the interaction between the two final 
state nucleons. Although it is perfectly true that the final 
state NN interaction significantly distorts the NN wave 
function but in the transition matrix element the con- 
tribution from all but the smallest NN separations will 
be strongly suppressed and the main effect may be at- 
tributed to the change of the normalization of the wave 
function at zero separation. If the non-interacting NN 
pair is described by a plane wave e"''" '', where k is the 
relative NN momentum {H = c = 1 units are used here- 
after), to account for final state interaction the latter 
must be replaced in the transition matrix element by the 
complete NN wave function (fc, satisfying outgoing 
spherical wave boundary condition at infinity. Neverthe- 
less, for a point-like interaction, we may set 

*~(fc,r)t w e-^'^ '' C(fc) (1) 



in the matrix element so that the final state interaction 
will be accounted for by multiplying the transition matrix 
element by the enhancement factor, defined as 

C(A:) = lim r)/e-* (2) 

The factor |C(fc)p that appears in the cross section rep- 
resents the ratio of two probabilities: one of finding 
the interacting NN pair at zero separation, while the 
other probability is associated with non-interacting par- 
ticles. By construction, when the final state interaction is 
turned off, the enhancement factor will be equal to unity. 
Expanding both, the numerator and the denominator on 
the right hand side of Q in partial waves, we have 

CW = lim^^-°(^^+^)-'^-^-(^'-)/-^!(^-^), (3) 

where ipe{k,r) ^ r^+^ for small r, ji{kr) is spherical 
Bessel function and Pe{k ■ r) denotes Legendre polyno- 
mial. Clearly, in the limit r — > in all higher partial 
waves will be suppressed by the centrifugal barrier, and 
only the contribution the from s-wave survives. Thus, we 
obtain a simple formula 

C{k)^Mk,Oy/k, (4) 

where prime denotes derivative with respect to r, and, as 
apparent from , the enhancement factor is determined 
by the slope of the wave function at the origin. To find 
this slope we must know the NN s-wave interaction and 
for simplicity we shall in the following assume that the 
latter takes the form of a spherically symmetric radial 
potential. The shape of this potential may be arbitrary 
but it must be of a short range. Given the NN poten- 
tial, we can integrate outward the appropriate wave equa- 
tion, containing both the nuclear and the Coulomb po- 
tential, generating numerically a regular solution uo(fc, r) 
(i.e. vanishing at the origin) whose derivative satisfies the 
boundary condition 

uo(fc,0)' = Co(?7)fc, (5) 

where i] denotes the Sommerfeld parameter and Co{rj)'^ — 
27r?7/[exp(27r7y) — 1] is the Coulomb barrier penetration 
factor. The sought for physical solution ipQ{k,r) occur- 
ring in 10} which is also regular, is necessarily propor- 
tional to uo{k,r), and, more explicitly, we have 

Mk,r)^[C{k)/Co{v)]uo{k,r). (6) 

Now, all we need to calculate C{k) is the asymptotic 
expression for the physical wave function. For r = R with 
R much bigger than the range of the nuclear potential, 
the physical wave function takes the form 

Mk-.R) = [C{k)/Co{v)]m{k,R) = 

^Fo{r,,kR) + fo{k)H+iri,kR), (7) 
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where H^{-q,kR) = GQ{r],kR) + iFo{rj,kR) with 
Go{ri,kR) and FQ{r],kR) being the standard Coulomb 
wave functions defined in and fo{k) = sin5e'* de- 
notes the s-wave scattering amphtude with S being the 
s-wave phase shift. The differentiation of |(7|) with respect 
to R, provides us with a second condition for the deriva- 
tives but it should be noted that uo{k,R) and u{k,Ry 
occurring in these two matching conditions are to be re- 
garded as known quantities. Indeed, they are fully speci- 
fied by the boundary condition at the origin ^ and can 
be either calculated analytically, or obtained by numer- 
ical methods. Therefore, what we end up with are two 
algebraic equations in which the two unknowns are the 
enhancement factor C{k) and the scattering amplitude 
fo{k) and the respective solutions, can be conveniently 
written as 



and 



C{k) 



w[H+{Tj,kR),uo{k,R)y 



w[F„{r],kR),uo{k,R)] 
w[H+irj,kR),uoik,R)Y 



(8) 



(9) 



where the symbol w[f,g] denotes the Wronskian defined 
as w[f,g] = fg' ~ f'g. We wish to recall that the spe- 
cific Wronskian involving the regular solution and the 
outgoing wave solution present in the denominator of 
and lO has been referred to as the Jost function 16]. It 
is worth noting that unlike the scattering amplitude © 
which can be expressed solely in terms of the logarithmic 
derivative of the regular solution, the enhancement factor 
depends separately upon both, u and u' , and there- 
fore the calculation of C(fc) requires off-shell information. 
In particular, a phase shift equivalent transformation of 
the potential may change the enhancement factor by or- 
ders of magnitude. Formula may be cast to the form 



C{k) = CwMik) 



{-k^)C,{vf 



w[Fo{v,kR),uo{k,R)y 



where 



CwAiik) 



' sin (5 



kCoiv) 



(10) 



(11) 



is the familiar Watson-Migdal factor pq|, depending 
solely upon the on-shell quantities, whereas the Wron- 
skian occurring in Hl()|l represents the off-shell correction. 

Formula ||SJ| deserves some further comments. It is easy 
to check that when both, the Coulomb and the strong 
interaction are switched off, the enhancement factor 
goes to unity. When the nuclear potential alone is set 
equal to zero, we have uo{k, R) = Fo{r], kR) and © yields 
C(fc) — Co{ri). Clearly, these are the right limits but they 
could not have been obtained with the Watson-Migdal 
factor alone which implies the importance of the off-shell 
correction. Nevertheless, close to threshold the off-shell 



factor varies slowly with energy, and the Watson-Migdal 
factor usually makes a good approximation. It should 
be kept in mind, however, that in result of multiplying 
the cross section by the Watson-Migdal factor the overall 
normalization is lost precluding an absolute calculation. 

Some authors choose to simplify further the Watson- 
Migdal enhancement factor (|ll|l by applying a Coulomb 
modified effective range formula for the phase shift. 
Thus, for the pp case, the phase shift can be obtained 
e.g. from the expression 



Cl^ri) fc cot(5 + 2kr]h{r)) 



-1/a + bk^/2 -Pk^/{1 + Qk^), 



(12) 



with h{r]) = Re ■0(1 + ^v) ^ log(^) where is the log- 
arithmic derivative of the gamma function. In H12() 
a = —7.83 fm and b = 2.8 fm denote, respectively, 
the experimental pp scattering length and the effective 
range and the remaining two parameters (P = 0.73 fm^ 
and Q = 3.35 fm^) are related to a specific NN poten- 
tial [221 . The approximation (|12l) has been popular in 
ceratin quarters and used in conjunction with Watson- 
Migdal formula (|ll|l for not very large k, results in a 
good approximation to the enhancement factor owing to 
a rather fortuitous cancellation of errors associated with 
different approximations. 



B. Examples 

We are going now to illustrate the results obtained 
in the preceding subsection by explicit calculations car- 
ried out for three phenomenological NN potentials which 
are: the delta-shell Vp , a Gaussian Vg , and the soft-core 
Reid Vr potential Among these three potentials 

only Reid potential has a repulsive short-range compo- 
nent while the two remaining ones are purely attractive. 
The reason for selecting these particular shapes is that 
they exhibit different behavior in the close to the ori- 
gin region: the delta-shell potential vanishes for small r, 
while a Gaussian potential shows maximum strength at 
r = 0, and, finally, the Reid potential is singular at the 
origin, i.e. Vr —^ oo. Our purpose is to examine how 
this very different off-shell behavior influences the en- 
hancement factor properties. Thus, as our first example 
we shall consider the delta-shell potential specified by a 
range R and a dimensionless depth parameter s, defined 
by the formula 



MVnir) = -{s/R)S{r-R), 



(13) 



where the values of the parameters are R = 1.84 fm 
and s = 0.906. Although this potential is not realistic, 
its main advantage lies in its simplicity: both the phase 
shift and the enhancement factor may be for this case 
obtained in an analytic form. A simple calculation gives 
the phase shift 



tan 6 



{s/p)Fo{v,P? 



l-(s/p)Fo(?7,p) Go(r/,p) 



(14) 




FIG. 1: Momentum dependence of So phase shift for differ- 
ent pp potentials. 



FIG. 2: Enhancement factor C(fc)p vs. k for different pp 
potentials. The dot-dashed curve has been obtained form 
itni and 



where p — kR, and, respectively, the enhancement factor 



l-(s/p)Fo(7;,p)i7o+(77,p)' 

It is apparent from IjlSfl that when the nuclear potential 
is switched off by setting s equal to zero the enhancement 
factor reduces to the Coulomb factor Co{r])- 

The Gaussian potential Vc{r) = Vq exp {—r'^/B?) has 
also two parameters, the depth Vq and the range i?, whose 
values are Vq = —31 MeV and i? = 1.8 fm and in this 
case the the solution of the wave equation will be ob- 
tained numerically. Finally, we consider the fully real- 
istic Reid soft core-potential [23|, given in the form of 
a superposition of three Yukawa potentials, which also 
requires numerical treatment. 

The results of our computations are presented in 
Figs, n and [21 In Fig. ^ we show the calculated pp 
phase shift vs cm. momentum as obtained from the 
different potentials. They are compared with the data 
and we can see that up to about 150 MeV all models 
agree well with experiment. For bigger momenta the sit- 
uation is less satisfactory, except for the Reid potential 
whose performance is still very good. Clearly, among the 
considered potentials only the Reid potential has a re- 
pulsive core and therefore is capable of reproducing the 
crossover from positve to negative values in the phase 
shift at 340 MeV. In Fig. [21 we compare the enhance- 
ment factors calculated from formula H1U|) for different 
potentials. We stipulated the normalization as to get in 
all cases the same height at the maximum. It is evident 
from Fig. [21 that very good agreement between differ- 
ent potential models results sustains for momenta up to 
about 150 MeV indicating that, apart from normaliza- 



tion, different off-shell properties have little impact on 
the shape of the enhancement factor. This can be easily 
understood since close to threshold effective range ap- 
proximation is usually adequate and all NN potentials 
are devised in such a way that they are capable of re- 
producing the experimentally determined effective range 
parameters. The main consequence of conducting a phase 
shift equivalent transformation is that the asymptotic 
wave function and its derivative acquire a common fac- 
tor preserving thereby the shape of the enhancement fac- 
tor. The differences appearing at larger momenta are not 
surprising as the different models do differ there also on- 
shell, as apparent from Fig.^ In particular, as seen from 
Fig. [21 the effective range formula makes a reliable ap- 
proximation for k not bigger than about 80 MeV. With 
the current high precision data, however, this approxi- 
mation is insufficient causing a too rapid fall-off of the 
enhancement factor. Indeed, already for the excitation 
energy as low as Q=15.5 MeV and the maximal momen- 
tum kmax = 120 MeV, \CwM{kmax)\^ Calculated from 
(|12|l is by a factor of 2 smaller from \C{kmax)\'^ obtained 
from (|10|l . with both functions normalized to yield the 
same height at the peak. In consequence, with the on- 
shell approximation the relative momentum distribution 
in the pp pprj reaction has a too rapid fall off. 



C. The cross section of the reaction pp pprj 

Since close to threshold only a small number partial 
waves contribute to the pp pprj transition amplitude, 
it is feasible to expand them in terms of angular momen- 
tum. The transition amplitude is an operator in spin 
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TABLE L List of allowed transitions ^"-+^1., '^"s+^lj^X 
for the lowest partial waves in the reaction pp — » pprj. 



even I 


oda I 


'Po ^^So,s 


^So 'Po,s 


'P2 -> ^So,d 


^D2 'P2,S 


'F2 ^So4 


'Po -> ^Pi,p 


'P2 ^D2,S 


'P2 'P2,P 


'F2 ^D2,S 


'F2^'Pl,p 




'F2~*'P2,P 




'Pi ^ 'Po,p 




'Pl^'Pl,p 




'Pl~''P2,P 




'F,^'P2,p 



partial wave amplitudes are functions of both, fc^ and 
but since these two quantities are linked by energy 
conservation one of them is redundant. Terms linear in 
k or q will be absent which can be seen as follows. The 
momentum dependence in a partial wave amplitude orig- 
inates from the spherical Bessel functions j t{kS,) j \{qif) 
occurring in the appropriate overlap integrals where ^ 
and r/ are the Jacobi coordinates canonically conjugated 
with k and q, respectively. For even orbital momentum 
the spherical Bessel function is an even function of the 
argument so that in the expansion of the partial wave 
amplitudes only even powers of the momenta will be ad- 
mitted. However, for practical reasons, it does not seem 
feasible to go beyond the second order in k and q. 

The above considerations lead us to take the following 
simple expression for the pp — > pp-q reaction cross section 



space that depends upon the three cm. momenta which 
determine the kinematics of the reaction, namely the ini- 
tial proton momentum p, the relative momentum of the 
final state protons fc, and, finally the momentum of the 
eta relative to the pp pair q (other possible choices will 
be discussed later on) . The quantum numbers associated 
with the initial state are: the angular momentum L, the 
total spin Si and the total angular momentum J. In the 
final state, the angular momentum and the the total spin 
of the pp pair are denoted as i and s / , respectively, and 
the relative motion of the eta is described by the angular 
momentum A. Among the above quantum numbers only 
J is conserved while the remaining quantum numbers are 
constrained by parity conservation and by the Pauli prin- 
ciple. Since the eta is a pseudoscalar meson, parity con- 
servation requires that L + i + \ must be an odd number. 
A list of possible transitions consistent with the above 
restrictions involving the lowest angular momenta is pre- 
sented in Table m These transitions can be classified as 
even or odd according to the value of the angular momen- 
tum £ of the final state protons. In the transition matrix 
element the even (odd) partial wave amplitude will be 
multiplied by an appropriate projection operator which 
is symmetric (antisymmetric) under the transformation 
fe — > — fe. Since the final state protons are indistinguish- 
able all observables must be invariant under the inter- 
change of the proton momenta, i.e. under the k —>■ k 
transformation. This means that all interference terms 
between even and odd states which are antisymmetric 
under fc — > — A; transformation are bound to vanish if one 
wants to respect Pauli principle. In consequence, inter- 
ference is allowed only between partial waves belonging 
to the same group which significantly reduces the number 
of terms in the cross section. 

For small excitation energy the final state pp inter- 
action appears to be the dominant effect and therefore 
it seems justified to neglect in the cross section all terms 
which do not contain the enhancement factor C(k). With 
this assumption the only contribution to the cross section 
comes from the even partial waves listed in Table |l| The 



da/dLips=\C{k)\^ [a + b P2{p ■ q)] + 

+Co{f]) [dReC{k) + elmC{k)]P2{p-k). (16) 

In H16|) we are using the standard notation where dLips 
stands for the invariant three-body phase space element. 
In (|16|l a denotes the modulus squared of the sole pro- 
duction amplitude which survives at threshold, associ- 
ated with the transition ^Pq — j- ^6*0,5, b represents the 
interference term between the latter amplitude and the 
•^Pq ^5*0, amplitude, and, the (d,e) pair describes, 
respectively, the interference with the '^P2 —>■ ^D2,s am- 
plitude. All of the above mentioned coefficients are func- 
tions of the final state momenta. Taking into account 
only threshold behavior, we can see that a is constant, 
b will be proportional to and (d,e) both are of the 
order of k^. In this situation, we have to expand a to 
the same order setting a = uq + ai q^ where ag and ai 
are two unknown parameters. The parameter oq can be 
absorbed in the normalization constant adjusted to the 
experimental total cross section value, and the ultimate 
expression for the cross section reads 

-^(x|C(fc)p + [x + yP2(p-9)]j + 

+ T^Co(r7) [zr ReC{k) + ImC(fc)] P2(p • fc), ^^^^ 

^max 



where x, y, Zr, Zi are real dimensionless parameters to be 
determined. For reasons of convenience we have scaled 
q and k by dividing them by their maximal values qmax 
and kmax, respectively. The parameter x represents the 
correction of the order q^ to the dominant transition, 
the parameters y and (z^, Zi) provide the measure of the 
admixture of the d- waves. 




FIG. 3: Spp distribution is compared with the data form 
Moskal et al. 0- The dotted curve corresponds to a constant 
^Po ^So s transition amphtude with on-sheU enhancement 
factor mil calculated from 1121 1. In the dashed curve the com- 
plete |C(fc)|^ has been used while in the full curve, at the top 
of that, a correction linear in energy has been included. 

III. COMPARISON WITH EXPERIMENT 

In Fig. |31 and in Fig. ^ we present the experimental 
distribution of the square of the effective pp mass Spp 
and the square of the rjp mass s^p, respectively. For a 
start, we assumed a constant production amplitude set- 
ting X — y = Zr = Zi = 0. Using the on-shell enhance- 
ment factor |Cvi/M(fc)P specified in (|ll|l with the effective 
range approximation (|12|) . from (|17|l we obtain the the- 
oretical distributions depicted by a dotted curve in Figs. 
13 and 21 Qualitatively, the dotted curves in Figs. El and 
0]are roughly in accord with experiment but they are far 
not satisfactory in quantitative terms. In both figures the 
peak in the calculated curves is too big and partially re- 
sponsible for that is the oversimplified enhancement fac- 
tor. As we already know, \CwM{k)\'^ calculated by using 
effective range formula (|12f) exhibits a too steep a fall 
and by normalizing the distributions to the total cross- 
section, to compensate that, the curves are moved up so 
that the peak gets magnified. Therefore, the simplest 
remedy is to abandon the on-shell approximation and 
from now on we will use the complete enhancement fac- 
tor calculated from (|10|l . The appropriate distributions 
are presented by the dashed curves in Figs. and ^ Al- 
though the resulting corrections go in the right direction 
bringing the calculation closer to experiment, but this is 
still not enough for providing a full understanding of the 
data. In this situation, it is interesting to know how much 
we can improve the agreement by disposing of the various 
corrections discussed in the preceding section and repre- 



FIG. 4: Srip distribution is compared with the data from 
Moskal et al. Q. The dotted curve has been obtained keeping 
only the dominant ^Po ^So s transition amplitude which is 
assumed to be a constant with |Cvi/M(fc)P calculated from 
11211 . In the dashed curve the complete |C(A:)p has been used 
while in the dot-dashed curve, at the top of that, a correction 
linear in energy has been included. The full curve accounts 
for also the interfernce with — > ^So d amplitude. 

sented by the adjustable parameters x, y, Zr, Zi. Before 
proceeding, however, it should be observed that in the 
case of the pp effective mass distribution a major simpli- 
fication takes place because, as apparent from (|17|l . af- 
ter integrating over the angles only the s-wave amplitude 
survives. Therefore, the parameters y, z^, Zi representing 
interference terms never enter Spp distribution. By con- 
trast, in the case of s^p distribution the interference terms 
in general do not vanish in result of angular integration. 
Since a non-relativistic approach is here completely ade- 
quate, in order to obtain the cross section as a function 
of Srip we have to introduce an equivalent Jacobi frame 
taking as the two independent momenta the 77-p relative 
momentum fci and the momentum of the other proton 
relative to this pair qi . The transformation between the 
two Jacobi frames, takes the form 

C) ^ {r :;^) CO 

where /i is the eta-p reduced mass and is the reduced 
mass of the eta and the pp pair. Clearly, substituting 
(|18|l in H17|l and integrating over the angles of fci and qi, 
the interference terms give non-vanishing contribution. 

Since the Spp distribution depends solely upon x, we 
adjusted this parameter by fitting the theoretical Spp dis- 
tribution to the experimental data with the best fit value 
of X being x=-0.514. The resulting distribution displayed 



in Fig|31by the full curve agrees now quite well with ex- 
periment. Actually, since the interference terms drop out 
and the 77P interaction gives here a small effect, there is 
really not much room for improvement other than includ- 
ing the correction in the the dominant ^Pq — )■ ^Sqs 
transition amplitude. Next, with the value of x in hand, 
we calculated the Sr,p distribution, still leaving out the 
interference terms i.e., setting y ~ Zr ~ zi = Q. The re- 
sulting cross section which does not involve adjustable 
parameters any more is presented in Fig01 by a dot- 
dashed curve. We can see that also in this case the 
agreement with experiment looks much better. The in- 
terference terms in (I17|l involving the parameters (z^, Zi) 
are probably much smaller than the term proportional to 
y because the former terms exhibit only a linear depen- 
dence upon the enhancement factor, so we ignored these 
terms adjusting the single parameter y to the experimen- 
tal Srfp distribution. With x fixed, the best fit value of 
y was y=3.38. The corresponding cross section is dis- 
played in Fig. 21 by the full curve and the agreement is 
quite good. We tried also to vary the parameters (zr, Zi) 
but this was not very successful as the fitting procedure 
attempts to reproduce a structure at the high energy end 
of the spectrum. With three adjustable parameters it is 
quite easy to produce a two-peaked distribution with a 
good . Since it is not quite certain that the data really 
reveal a two-peaked distribution we did not pursue this 
fit any further. It should be also mentioned that when 
the interference term linear in C(k) is accounted for, the 
calculation becomes very sensitive to the NN potential 
used to obtain C(k) because the normalization of C(k) 
can not be absorbed in the overall normalization of the 
cross-section. 

In Fig. [5I we are going to present the changes in the 
excitation energy distribution caused by adding an en- 
ergy dependent term to the eta production amplitude 
(the interference terms in (|17|l give no contribution) . The 
case x=0 corresponds to a constant production ampli- 
tude (dotted curve in Fig. El, while x=-l (dashed curve 
in Fig. [SJ represents the lower limit for this parameter 
which must be imposed to prevent the cross section from 
going negative. The best fit value x«-0.5 necessary to re- 
produce the effective mass distributions for Q=:15.5 MeV 
is roughly midway between these values and the corre- 
sponding cross section is presented by the full curve in 
Fig. El The linear energy dependence introduced in the 
amplitude of the eta production has little impact on the 
total cross section. In fact, the two curves where x<0 ex- 
hibit slightly better agreement with experiment as com- 
pared with the dotted curve corresponding to constant 
production amplitude. To obtain our curves we used the 
full expression for the enhancement factor. We wish to 
note that as a result of using the on-shell approximation 
for the enhancement factor the calculated cross section 
is underestimated at large Q. 

With a good fit to the effective mass distributions for 
Q=15.5 MeV, it is interesting to know whether the ap- 
proach presented above makes sense for Q=41 MeV since 
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FIG. 5: The total cross section for the pp pprj reaction 
versus the excitation energy Q for different values of the pa- 
rameter X. The experimental data are from 0]. 



for this value of the excitation energy the experimental 
data are also available in the literature . Strictly speak- 
ing, for this much larger excitation energy the validity of 
the simple formula p7|l becomes questionable and the in- 
clusion of p-wavcs might be indispensable. Nevertheless, 
it is useful to find out what a simple model can do for 
larger Q. The results of our computations are presented 
in Figs. El and where they are compared with the data 
from 5] . In the calculation we neglected the interference 
terms, setting ?; = = = in (|17|) . The dotted curve 
in Fig. Eland Fig.[3corresponds to a constant eta produc- 
tion amplitude. Whilst the Srip distribution is not too far 
off the experiment, the Spp distribution is in a bad shape 
as the peak caused by pp FSI is definitely much too big. 
The dashed curve is obtained by adopting for x the same 
value that for Q=15.5 MeV gave the best agreement with 
experiment i.e., we are pretending that this parameter 
does not vary with Q. This brings improvement in both 
cross sections and in fact the s^p distribution agrees with 
the data quite well. Since the parameter x may depend 
upon Q, and for Q=4I MeV its value could be different 
than for Q=I5.5 MeV, we allowed x to vary, adjusting its 
value to the data from This time the best fit value 
turns out to be x=-0.863 and the corresponding distri- 
butions are presented by full curves in Figs. El and It 
is not possible to improve the agreement in both distri- 
butions using the same value of the parameter x, and 
therefore the fit depicted by the full curve is a compro- 
mise. Although the corrections go in the right direction 
but the agreement is still not satisfactory and we have 
to accept the fact that for Q=41 MeV the assumed func- 
tional form of the cross-section is incomplete. We tried 
to include the interference terms with the d-waves but 
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the fit was unsuccessful. Clearly, further extensions call 
for additional parameters but this does not seem to be 
affordable with the present data, and therefore we stop 
here. 

Summarizing, a simple model in which the pp — > pprj 
cross section contains only the dominant •^Pq ^So,s 
partial wave amplitude corrected by the final state pp 
interaction is capable of explaining the current experi- 
mental data at an excitation energy Q=15.5 MeV. The 
experimental pp and the rjp effective mass distributions 
can be both reproduced if: (i) in the eta production am- 
plitude a correction linear in energy is admitted and, (ii) 
the full pp FSI enhancement factor without the on-shell 
approximation is used. In the presented calculation the 
Typ final state interaction has not been included explicitly 
for reasons outlined in the Introduction but since the pa- 
rameters X and y have been adjusted to the data, they 
effectively account for this effect. 



FIG. 6: Spp distribution for an excitation energy Q=41 MeV 
is compared with the data form Abdel-Bary et al. 1^. In 
our calculation we have neglected all higher waves except the 
^Po ^So s transition amplitude. The energy dependence in 
the latter is proportional to the parameter x, cf. 1171 . 
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FIG. 7: Srjp distribution for an excitation energy Q=41 MeV 
is compared with the data form Abdel-Bary et al. |3- In 
our calculation we have neglected all higher waves except the 
^Po ^So s transition amplitude. The energy dependence in 
the latter is proportional to the parameter x, cf. 1171 . 



